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ABSTRACT 



We present the results of observations of the radio emission from Supernova 
2001gd in NGC5033 from 2002 February 8 through 2006 September 25. The data 
were obtained using the Very Large Array at wavelengths of 1.3 cm (22.4 GHz), 
2 cm (14.9 GHz), 3.6 cm (8.4 GHz), 6 cm (4.9 GHz), and 20 cm (1.5 GHz), with 
one upper limit at 90 cm (0.3 GHz). In addition, one detection has been provided 
by the Giant Metrewave Radio Telescope at 21 cm (1.4 GHz). SN 2001gd was 
discovered in the optical well past maximum light, so that it was not possible 
to obtain many of the early radio "turn-on" measurements which are important 
for estimating the local circumstellar medium (CSM) properties. Only at 20 cm 
were turn-on data available. However, our analysis and fitting of the radio light 
curves, and the assumption that the Type lib SN 2001gd resembles the much 
better studied Type lib SN 1993J, enables us to describe the radio evolution as 
being very regular through day ~550 and consistent with a nonthermal-emitting 
model with a thermal absorbing CSM. The presence of synchrotron-self absorp- 
tion (SSA) at early times is implied by the data, but determination of the exact 
relationship between the SSA component from the emitting region and the free- 
free absorption component from the CSM is not possible as there are insufficient 
early measurements to distinguish between models. After day ~ 550, the radio 
emission exhibits a dramatically steeper decline rate which, assuming similarity 
to SN 1993J, can be described as an exponential decrease with an e-folding time 
of 500 days. We interpret this abrupt change in the radio flux density decline 
rate as implying a transition of the shock front into a more tenuous region of 
circumstellar material. A similar change in radio evolution has been seen earlier 
in other SNe such as SN 1988Z, SN 1980K, and SN 1993J. 

Subject headings: supernovae, individual (SN 2001gd), galaxies, individual (NGC5033), 
radio continuum, stellar evolution: massive stars 



Introduction 



Su pernova (SN) 2 001gd was discovered at magnitude 16.5 on 2001 November 24.820 
UT by iNakano et al\ ( I200ll ). located ~ 3' north-northwest of the center of NGC 5033. 
Nakano et all toOlh report a position of a = 13 ft 13 m 23!89, 5 = +36 38'17?7 (J2000), for the 
SN and note that nothing was detected at this position in previous CCD i mages taken 179 
days earlier (2001 May 29) to a limiting magnitude of 19. Ten days later, iMatheson et al. 
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( 120011 ) obtained a spectrum of the SN, identifying it as a Type lib SN well past its maximum 
light, with weak helium lines and nebular-phase lines of magnesium, oxygen, and calcium 
beginning to dominate. They noted that this spectrum was nearly identical to that of Type 
lib SN 1993J, on day 93 after its explosion. Using this suggestion to constrain the explosion 
date for the SN, we assume that the optical discovery occured ~ 83 days after the explosion 
and, therefore, for calculating the age of the supernova we as sume an explosion date of 2001 
September 03, which is the same explosion date assumed by iPerez- Torres et all ( 120051 ). 



Tullvi (119881) gives the distance to NGC 5033, as 18.7 Mpc (calculated with H 



75 km s _1 Mpc" 1 ). Accepting their distance calculation and correcting it with an updated 
value of the Hubble constant H = 65 km s _1 Mpc -1 gives a distance to the SN of ~ 21.6 Mpc. 
We note, however, that NG C 5033 is the host of several historic al supernovae, includin g the 
radio emitti ng SN 1985L jVan Dvk et a,l\\l99^ }. If we follow Ivan Dvk etaH jl99Rh and 
estimate the iTullyi ( 19881 ) distance with Hubble constant Hq = 70 km s _1 Mpc" 1 , we obtain 
a distance to NGC 5033 of ~ 20 Mpc. Since a change in the distance estimate of < 10% 
will not significantly affect any of our estimates of the SN physical properties, we will, for 
consistency with the previous ra dio work, use in this paper the Van Dyk distance of 20 Mpc. 
However, it should be noted that IPerez- Torres et all (120051 ) assumed a much shorter distance 
of 13.1 Mpc. 



SN 2001gd was first detected at radio wavelengths on 2002 February 8.54 by lStockdale et al. 

( b()02h . using the Very Large Array (VLA)Q at 1.3 cm (22.4 GHz), 3.6 cm (8.4 GHz), 6 cm 
(4.9 GHz), and 20 cm (1.5 GHz). They report a radio position of a = 13 /l 13 m 23!889, 
5 = +36°38'18 / . / 14 (J2000) with uncertainties o f ±0'/ 2 in each coordinate, which agrees with 
the reported optical position of iNakano et al\ (120011 ) to within the errors. 



Stockdale et all ( 120031 ) presented the results of initial radio monitoring of SN 2001gd 



between 2002 February 8 and 2002 October 28. Their results showed that the radio emission 
evolves very smoothly, consisten t with the nonthermal-emitting, thermal-absorbing model 
described by IWeiler et al\ ( 120021 ). Their results also indicate a CSM with significant struc- 
ture embedded in a uniform medium, similar to that proposed for Type lib SN 1993 J 
( IVan Dyk et aZ.I I1994J : IWeiler et all 120071 ). The similar optical spectra and early radio evo- 
lution of both SN 1993J and SN 2001gd indicate that the two Type lib SNe may have a 
number of similarities. 



Perez- Torres et all (120051 ) conducted VLBI observations of SN 2001gd at 3.6 cm (8.4 



GHz) on 2002 June 26 and 2003 April 8. Although they were unable to resolve any structure 



1 The VLA telescope of the National Radio Astronomy Observatory is operated by Associated Universities, 
Inc. under a cooperative agreement with the National Science Foundation 
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in the SN, they were able to estimate the SN angular size through interferometer visibility 
fitting and found that their data suggested a possible deceleration in the expansion (although 
their data were not conclusive). Using o ur assumed distance of 20 Mpc to the host galaxy, 



and the angular size limits estimated by iPerez- Torres et al\ (120051 ). the upper limit to the 



expansion velocity of the blastwave is < 22, 000 km s _1 for the 2002 June 26 VLBI observation 
(day 296). Depending on the exact nature of the emitting region, the upper limit may be 
as low as < 19,000 km s _1 . Thus, for the purposes of this paper, and for comparison with 
SN 1993 J, we assume an expansion velocity of 15,000 km s _1 . 



Perez- Torres et al\ (120051 ) conducted multifrequency radio observations with the VLA 
on 2003 April 8 to complement their VLBI data, obtaining a spectrum consistent with 
synchrotron emission partially suppressed by free-free absorption. They also obtained X-ray 
data on SN 2001gd from XMM-Newton observations of NGC5033 on 2001 July 2 and 2002 
December 18 and report the detection of X-ray emission on the latter date consistent with 
the expected reverse shock emission from typical Type II SNe. 

Recent VLA observations indicate a dramatic change in the evolution of the radio emis- 
sion from SN 200 lgd. These new observations and their implications are presented in this 
paper. 



2. Radio Observations 



New radio observations of SN 2001gd have been made with the VLA at 1.3 cm (22.4 
GHz), 2 cm (14.9 GHz), 3.6 cm (8.4 GHz), 6 cm (4.9 GHz), and 20 cm (1.5 GHz), with one 
upper limit at 90 cm (0.3 GHz) from 2002 November 08 through 2006 September 25 and 
are presented in Table [TJ and Figure H For completeness a nd ea se of reference, we i nclud e 
the previously published results from lStockdale et all (120031) . and IPerez- Torres et all ( 120051 ). 



as well as one measurement from the GMRT (IChandra et al. 1 120021 ) , in both Table [Tj and 
Figure [TJ 

The techniques of observation, editing, calibration, and error estimation are described 



in previous publications on the radio emission from SNe (see, e.g., IWeiler et q/.lll986l . 12002 



Sramek fc Weiler 1 120031 ) . We used 3C286 as the primary calibrator for flux density boot- 
strapping, while J1310+323, with a defined position of a = 13 h 10 m 28!66, 5 = +32°20 / 43'. / 8 
served as the phase calibrator and secondary flux density calibrator. The flux densities of 
J1310+323 are shown in Table [2] and plotted in Figure [2j Like most phase calibrators, 
J1310+323 shows some variation in its flux density because phase calibrators are selected for 
their proximity to the SN and unresolvable angular size rather than their flux stability. The 
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primary calibrator is therefore used to determine the absolute flux density of the secondary 
calibrator at each epoch of the observations. 

The listed flux density measurement errors are determined based on both the intrinsic 
RMS background noise level of the particular radio map, and a fractional error, e, propor- 
tional to the SN flux. The RMS map noise results from both small unresolved fluctuations 
in the background emission and random map fluctuations due to receiver noise, while the 
fractional error, e , accounts for the inaccuracy of VLA flux density calibration (see, e.g., 



Weiler et a/.lll986l ) and possible deviations of the primary calibrator from an absolute flux 
density scale. The final errors (cr/) given for the measurements of SN 2001gd are taken as 
the addition in quadrature of these components: 



a) = (eSof + *l (1) 

where Sq is the measured flux density, o"o is the map RMS background noise for each obser- 
vation, and e = 0.15 for 90 cm, e = 0.1 for 20 cm, e = 0.05 for 6 cm and 3.6 cm, e = 0.075 
for 2 cm, and e = 0.1 for 1.3 cm. 



3. Parameterized Radio Light Curves 



From the most recent RSN modeling discussion of lWeiler et al\ (120021 ) and lSramek fc Weiler 



( 120031 ) . we use a parameterized model : 



V \ Q / t tn\ ( "\- 6 TCSM dumps \ / \ g internal 



5 ( m Jy) = *i tt^t U^M e " Tcxtcrnal — — — ( 2 ) 

V5 GRzJ VI day; V ^csM clumps / V Wnai ' 



with 



^external — r CSM homogcncous + ^distant; (3) 

where 



v \- 2 - 1 ft-t xS 



TcSMho — K2 \5GRz) U day/ ' lj 

-2.1 

distant = K A ( ) (5) 



- 6- 



and 



-2-i ft-t 
5 GHz/ \ \ day 



r csM clumps = K 3 ( F7TEJI J ( TTJ:3 ) (6) 



with i^i, K 2 , K 3 , and determined from fits to the data and corresponding to the flux 
density (Ki), homogeneous (K 2 , K4), and clumpy or filamentary (K 3 ) free-free absorption 
(FFA) at 5 GHz, formally but not necessarily physically, one day after the explosion date 
to- For a more detailed explanation of the physical meaning of the parameters, the reader 
should consult the references cited above. 

Since it is physically realistic and may be needed in some RSNe where radio observations 
have been obtained at early times and high frequencies, Equation (j2J) also includes the 
possibility for an internal absorption term. This internal absorption (r interna i) term may 
consist of two parts - synchrotron self-absorption (SSA; r internalssA ) , and mixed, thermal 
FFA/non-thermal emission (r in t erna i FFA ). 

^internal 'HnternalggA 'HnternalppA V ) 



TinternalgsA -^5 



.. X' 1 



5 GHz/ V. 1 da y 



/ v \- 21 /t -t y 

rinternal FFA - ^6 {^^) [jd^ ) ^ 

with corresponding to the internal, non-thermal (v a ~ 2 5 ) SSA and Kq corresponding to 
the internal thermal (v~ 2A ) free-free absorption mixed with non-thermal emission, at 5 GHz, 
formally but not necessarily physically, one day after the explosion date t . The parameters 
5" and 5"' describe the time dependence of the optical depths for the SSA and FFA internal 
absorption components, respectively. 

Application of this basic parameterization has been shown to be effective in describ- 
ing the physical characteristics of the presupernova system, its CSM, and its final stages 
of evolution before explosion for objec ts ranging from the two decades of monitoring the 



complex radio emission from SN 1979C (jMontes et a/.ll2000l) through the unusual S N 1998bw 



(GRB980425) flWeiler. Panada. & Montes Il2001[ ) to the very well studied SN 1993J flWeiler et al. 



20071 1 . 
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4. Radio Light Curve Description 



The multi-frequency radio data for SN 200 lgd are listed in Table [T] and displayed in 
Figure [TJ Examination of the data indicates that radio evolution of the SN is best described 
by two distinct time intervals: an "early" interval before day ~ 550 and a "late" interval 
thereafter. Although the precise date of the observed "break" i s not w ell defined, it appears 
that around the time of the observations of iPerez- Torres et all (120051 ) on 2003 April 8 (day 
582) the rate of flux decline markedly steepened such that the data can no longer be ade- 
quately described by the model appropriate for earlier times. Simil ar "breaks" haye bee n 
observed in the radio light curv es of other RSNe , including SN 198 8Z ([Williams et a/.ll2002l ). 
SN 1980K jMontes et ai.lll998h . and SN 1993J (jweiler et alhoQ^ . 



An extensive search for the best parameter fits to equations 



,2 \ 



was carried out 

on the "early" interval (day < 550) by minimizing the reduced \ 2 (Xted)- Since there are 
only four radio detections after the measurements on day 582, the "late" radio light curve 
behavior is very poorly defined except to note that the few detections and the several upper 
limits indic ate that the declin e rate has significantly steepened and, following the example of 
SN 1993 J (IWeiler et all 120071 ). are consistent with an exponential decline with an e-folding 
time of 500 days. 

The early epoch fitting parameters are consistent with a homogeneous {K 2 ) CSM, al- 
though more complex structures such a CSM containing filamentary (-K3) components cannot 
be ruled out. No evidence is seen for any distant, homogeneous thermal absorption (K± = 0), 
and the very early data (< 100 days) are too sparse to determine from the fitting if SSA (K 5 ), 
or mixed, internal, free-free absorption/nonthermal emission (Kq) may have been present. 

However, we have estimated the brightness temperatures for the "early" (before day 
~ 550) 20 cm emission based on a pure FFA model and found tha t they exceed the limiting 
temperature for synchrotron em ission of ~ 10 1L5 K calculated by iReadhead I (11994J ). which 



was also the case for SN 1993 J (IWeiler et all 120071 ). Given our inability to model this early 
absorption based purely on our sparse "early" data, we have chosen to assume that SN 2001gd 
has a similar brightness temperature evolution (~ 10 11 K at peak) for the 20 cm emission as 
it becomes optically thin. Given that SN 2001gd and SN 1993J have similar characteristics at 
"early" times, we estimate, as was the case for SN 1993J, that a significant SSA component 
is required for SN 200 lgd along with the FFA component determined from fitting which 
describes the rapid rise of the 20 cm data well. The best parameters which we can estimate 
for the "early" epoch of SN 2001gd with these assumptions are given in Table [31 We note an 
apparent similarity between the 5 and 5" terms from our modeling of the SN 2001gd radio 
emission. This is similar to the case of SN 1993 J where the two values were 5 = —1.88 and 
5" = —2.05. However, it must be stated that, from pure modeling, our sparse "early" data 
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only constrain the values of 5 and 5" for SN 2001gd to lie between —1 and —2. To narrow 
the range of these values, we ha ve therefore assume d a value of 5 = —1.88 to match the value 
for 5 determined for SN 1993J JvVeiler et AoOt! ). 



The "late" period appears to be optically thin at all observed frequencies with a much 
steeper decline rate than is possible with the "early" data determined (3 = —0.92 (shown 
as the dashed line i n Figure ED). Again , if we assume a similarity to the exponential decline 
seen for SN 1993J (jWeiler et all 120071 ). the late epoch of SN 2001gd is consistent with an 
exponential decline after day 550, with an e-folding time of 500 days. This observed "break" 
in the radio light curves of SN 2001gd is interpreted as indicative of a rapid transition of the 
blastwave into a much more tenuous medium. This oc curs much ea r lier fo r SN 2001gd (day 
~ 550) than that seen for SN 1993J (day ~ 3, 100) bv IWeiler et ~aH J2007h . 



Before the break, SN 2001gd is still somewhat optically thick at 20 cm, but after the 
break it appears to be optically thin. This discounts the possibility of the break being the 
result of a sudden change in the shock velocity or magnetic field properties at the shock 
front since such changes would not yield an optical depth reduction. Given the similar 



expan sion velocities of < 19, 000 km s 1 for SN 2001gd (derived fro m (IPerez- Torres et al. 



2005) ; see the Introdu ction) and 11, 000 - 15, 000 km s" 1 for SN 1993J flMarcaide et al 



1997 



Marcaide et a/.ll2007l ). it is clear that both supernovae had a period of enhanced mass- loss 
rate preceding their explosions. Thus, assuming a blastwave velocity for SN 2001gd of 15, 000 
km s _1 and a "standard" pre-explosion wind velocity of 10 km s _1 , this enhanced mass-loss 
persisted for the ~ 1, 800 years prior to the explosion. However in the case of SN 1993J, 
this phase began ~ 8,800 years earlier than the explosion of its progenitor star. The shorter 
duration of the high presupernova mass- loss rate for SN 200 lgd is also borne out in its 
faster rise and higher luminosity, indicating there was less intervening, presupernova wind 
established, CSM to contribute to the FFA at early epochs. 



4.1. The Mass-Loss Rate for the SN Progenitor 

Through use of this modeling, a number of physical properties of SNe can be determined 
from the radio observations. One of these is the mass-loss rate from the SN progenitor prior 
to explosion. From the Chevalier (1982a,b) model, the turn-on of the radio emission for RSNe 
provides a m e asure of the presupernova mass-loss rate to wind velocity ratio (M/w win d) and 



Weiler et al\ (119861 ) derived this ratio for the case of pure, thermal, external absorption by 



a homogeneous medium. 



Given the limited amount of early radio data available for SN 200 lgd, there is no need to 
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invoke more complex models than the simple, uniform, external th ermal absorption model, 
so that r{t = 1 day) = K 2 and equation (16) of IWeiler et al\ (119861 ) can be written: 



M(M Q yr- 1 ) 
(^wind/10 km s" 1 ) 



3.0 x 10~ 6 T^Lr m" 

1 v - JL - >lvl homogcncous 



1.5 



ti 



45 days 



Vi 



10 4 km s" 1 



1.5 



1.5 



1.5m 



T 



10 4 K 



0.68 



(10) 



Here, the extra factor is a small correction that takes into account the fact that, for 
SN 2001gd, the CSM density behaves like p oc r~ 161 (see below) instead of r~ 2 as it would 
under the usual constant mass-loss rate assumption. The factor <fi is given by the square root 
of the ratio of the integration constant for r in the case of p oc r -L61 to the one appropriate 
for p oc r~ 2 , i.e. 



2 x 1.61 - 1 
2x2-1 



0.5 



0.86 



Since SN 2001gd is a Type lib supernova as SN 1993J was, and appears to be quite 
similar to it, we will adopt for SN 2001gd some of the same parameters as measured for SN 
1993J. Thus, for SN 2001gd we assume v\ = 15, 000 km s _1 at t = 45 days, which is a value 



consistent with the results of iPerez- Torres et all ( 120051 ) as discussed earlier 



We also adopt values of T = 20, 000 K, = 10 km s" 1 (which is appropriate for a 

RSG wind), t = (iecm peak — to) days from our bes t fits to the radio data, and m = 0.845, as 
measured for SN 1993 J by IMarcaide et all (120071 ). With the assumptions for the blastwave 
and CSM properties discussed above, and the results for the best-fit parameters listed in 
Table [3], our estimated presupernova mass-loss rate is M = 5.0 x 10~ 6 M yr^ 1 at 8.6 years 
before explosion. 

The fitting parameter value of 5 = —1.88 indicates that the CSM density behaves like 
p oc r _L61 so that the mass-loss rate for SN 2001gd was not constant but was decreasing in 
the years leading up to the explosion, i.e. M oc r 2 pw7 win d oc r 39 for a constant w w ind- 



15, 000 km s" 1 at t 



45 



Since the shock speed for SN 200 lgd is assumed to be 
days, but variable with time (v oc t" 0155 for m = 0.845 from IMarcaide et all 120071 ). and 
the pre-supernova wind velocity for an RSG is typically ~ 10 km s _1 , we calculate that the 
abrupt decline in the radio light curves starting around day ~ 550 implies this change took 
place in the presupernova stellar wind ~ 1,800 years before explosion when the mass- loss 
rate was as high as 3.9 x 10~ 5 M Q yr -1 . Integrating the mass-loss rate over the last ~ 1, 800 
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years before explosion, we find that, during that time, the progenitor star shed ~ 0.05 M Q 
in a massive stellar wind. At earlier epochs of the progenitor's evolution the mass-loss rate 
was considerably lower, as indicated by the radio light curve "break" discussed above and 
the transition of the shock wave to a lower density CSM at that time. 



Conclusions 



Figure 1 and Table 3 show that the radio ligh t curves for SN 2001gd can be described 
by standard RSN models jWeiler et alWimd . Il99pl . I2OO2I : iMontes. Weiler. fc Panagk]ll997l ) 
during its "early" evolution up to day ~ 550. Aft er day ~ 550 the f lux density decline rate 
clearly steepened to, in comparison with SN 1993 J JWeiler et oJ2007h . a form consistent with 
an exponential decay with e-folding time of 500 days. The "break" in the radio light curves for 
SN 2001gd is interpreted as due to a chan ge in the average CS M density. S imilar "breaks 



have been identified earlier in SN 1980 K (IMontes et aZ.lll998l ). SN 1988Z flWilliams et al. 
2002h . and SN 1993 J JWeiler et aZj|2007h . 



Normal assumptions for presupernova stellar wind properties imply a dense wind with 
mass-loss rate decreasing for the last ~ 1, 800 years before explosion from M ~ 3.9 x 10 -5 M & 
yr _1 to M ~ 5.0 x 10~ 6 M yr _1 around the time of explosion. The mass-loss rate was sharply 
lower for times earlier than ~ 1, 800 years before explosion. 



KWW wishes to thank the Office of Naval Research (ONR) for the 6.1 funding sup- 
porting his research. CJS is a Cottrell Scholar of Research Corporation and work on this 
project has been supported by the NASA Wisconsin Space Grant Consortium. NP is As- 
tronomer Emeritus at the Space Telescope Science Institute (STScI) that kindly provided 
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Table 1. Flux Density Measurements for SN 2001gd a 



Obs. Date 


Age 


Tel. 


S (90 cm) 


S (20 cm) 


S (6 cm) 


S (3.6 cm) 


S (2 cm) 


S (1.3 cm) 


CUT) 


(days) 




(mjy) 


(mjy) 


(mjy) 


(mjy) 


(mjy) 


(mjy) 


2001 Sep 03 


EE 
















2002 Feb 08 


159 


VLA-A 




0.85 ±0.15 


5.05 ± 0.26 


4.09 ±0.21 




1.37 ±0.18 


2002 Mar 02 


180 


VLA-A 




1.41 ±0.20 


5.41 ± 0.28 


3.73 ±0.20 


2.18 ±0.24 


1.34 ±0.18 


2002 Mar 14 


192 


VLA-A 




1.84 ±0.37 


4.46 ± 0.24 








2002 Mar 23 


201 


VLA-A 




1.83 ±0.36 




3.66 ±0.21 


2.33 ±0.32 


1.00 ±0.19 


2002 Apr 06 


215 


VLA-A 




1.99 ±0.37 


5.33 ± 0.28 


3.78 ±0.21 


2.66 ±0.30 




2002 Jun 12 


282 


VLA-B 




3.23 ±0.37 


4.84 ±0.25 


3.10 ±0.17 






2002 Jun 26 b 


296 


VLBA 








3.83 ±.19 






2002 Jul 16 


317 


VLA-B 




2.73 ±0.31 


4.26 ± 0.23 


2.39 ±0.14 






2002 Sep 22 c 


384 


GMRT 




3.55 ±0.22 










2002 Oct 28 


421 


VLA-C 








2.07 ±0.14 


0.89 ±0.18 


0.81 ±0.16 


2002 Nov 08 


432 


VLA-C 




4.22 ± 0.48 


3.26 ± 0.20 








2003 Apr 08 b 


582 


VLA-D 




3.89 ± 0.40 


1.85 ± 0.24 


1.02 ±0.05 


0.63 ±0.13 


0.39 ±0.12 


2003 May 23 


627 


VLA-A 


< 3.09 






0.96 ±0.10 




< 0.39 


2003 May 26 


630 


VLA-A 




3.47 ± 0.45 


1.67 ±0.10 








2004 Sep 10 


1104 


VLA-A 




1.04 ±0.14 




< 0.18 


< 0.79 


< 0.79 


2005 Jun 14 


1380 


VLA-C 




< 1.12 


< 0.21 


< 0.14 


< 0.48 


< 0.42 


2006 Jan 24 


1605 


VLA-D 






< 0.39 


< 0.6 






2006 May 29 


1729 


VLA-BnA 




< 0.32 










2006 Sep 25 


1849 


VLA-B 






< 0.29 









a All upper limits are 3cr. 

b From IPerez- Torres et~al Perez- Torres et al ] l|2005h also measured a 3<r upper limit of < 0.3 mjy for SN 200 lgd 

on 2003 Apr 08 at 0.7 cm (43.3 GHz). 

^Chandra et al. I d2002h 
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Table 2. Flux Density Measurements for Secondary Calibrator J1310+323 



Obs. Date 


Age 


Tel. 


S (90 cm) 


S (20 cm) 


S (6 cm) 


S (3.6 cm) 


S (2 cm) 


S (1.3 cm) 


CUT) 


(days) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 




2002 Feb 08 


159 


VLA-A 




1.70 


1.66 


1.77 




2.15 


2002 Mar 02 


180 


VLA-A 




1.70 a 


1.69 


1.75 


2.03 


2.16 


2002 Mar 14 


192 


VLA-A 




1.70 a 


1.70 a 








2002 Mar 23 


201 


VLA-A 




1.68 




1.80 a 


2.08 


2.32 


2002 Apr 06 


215 


VLA-A 




1.69 


1.70 a 


1.83 a 


2.07 




2002 Jun 12 


282 


VLA-B 




1.67 


1.73 


2.01 






2002 Jul 16 


296 


VLA-B 




1.72 


1.77 


2.05 






2002 Oct 28 


421 


VLA-C 








2.13 


2.55 


2.73 


2002 Nov 08 


432 


VLA-C 




1.70 


1.83 








2003 Apr 08 b 


582 


VLA-D 




1.54 


1.66 




2.40 


2.91 


2003 May 23 


627 


VLA-A 


2.74 






2.37 




2.93 


2003 May 26 


630 


VLA-A 




1.42 


1.65 








2004 Sep 10 


1104 


VLA-A 




1.14 




2.47 


2.35 


2.24 


2005 Jun 14 


1380 


VLA-C 




1.45 


1.81 


1.89 


1.76 


1.71 


2006 Jan 24 


1605 


VLA-D 






1.44 


1.35 






2006 May 29 


1729 


VLA-BnA 




1.58 










2006 Sep 25 


1849 


VLA-B 






2.42 









a Value was not measured but was interpolated from the closest calibration at that frequency. 

b From Ipirez-Torres et al\ J2005h ; IPerez-Torres et al\ J2005I) also measured a value of 2.58 Jy for J1310+323 on 
2003 Apr 08 at 0.7 cm (43.3 GHz). 
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Table 3. Fitting Parameters for SN 2001gd 



Parameter 


Early (day < 550) 




7.5 x 10 2 


a 


-0.94 


P 


-0.92 


K 2 


1.5 x 10 3 


5 


-1.88 


K 5 


1.7 x 10 2 


5' 


-1.5 




2001 Sep 03 


S6cm peak ( m Jy) 


7.96 


L 6c m peak (erg s _1 Hz" 1 ) 


3.8 x 10 27 


t = (U cm peak ~ *o) (days) 


80 


M (M jY- v f 


(0.5-3.9) x 10~ 5 


Assumed distance (Mpc) 


20 



a Defined in the Introduction. 

b Assuming w wind = 10 km s" 1 , v { = tWtwave = 
15, 000 km s _1 , T = 20, 000 K, and t ; = 45 days. 
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Fig. 1 — Radio "light curves" for SN 2001gd in NGC 5033. The five wavelengths, 1.3 cm 
(22.4 GHz; top left), 2 cm (14.9 GHz; middle left), 3.6 cm (8.4 GHz; middle right), 6 cm 
(4.9 GHz; bottom left), and 20 cm (1.5 GHz; bottom right) are shown together with their 
best-fit model light curves. The SN age is in days since the adopted explosion date of 2001 
September 03. Because the decline i ndex of the radio emission steepened after day ~ 550, 
following the example of SN 1993 J (jWeiler et all 120071 ). we have fitted the "late" data after 
day 550 with an exponential decline with an e-folding time of 500 days. The inverted open 
triangles represent 3cr upper limits and the dotted lines show the extrapolation of the fit to 
the "early" data before day ~ 550 given in Table [3] 



calibrator J1310+323 



5 F 

2 
1 
5 

2 
1 
5 

2 
1 
5 

2 
1 
5 

2 
1 



— 

- 1.3 cm 
L 1 1 1 1 1 1 


■ ■■ 


■ 

■ 


2 cm 
1 1 1 1 1 1 1 


1 1 — 1 1 1 1 1 1 

■ ■ 

■ ■■ 


1 Hz 

■ 

■ 


3.6 cm 
1 1 1 1 1 1 1 


1 1 — 1 1 1 1 1 1 

" " 


1 Hz 

■ 

■ 

■ 


6 cm 
1 1 1 1 1 1 1 


1 1 — 1 1 1 1 1 1 

■ ■■■ ■■ * ■■ 


1 Hz 

■ 

■ 

■ 


20 cm 
E 


1 1 — 1 1 1 1 1 1 

■ mmm ■ ■ ■ 

■ 




1 Hz 

■ ■ 

■ 



100 1000 
Time since explosion [days] 



Fig. 2. — Flux density measurements of the secondary calibrator J1310+323 at the five 
wavelengths of 1.3 cm (22.4 GHz; top), 2 cm (14.9 GHz; one below top), 3.6 cm (8.4 GHz; 
middle), 6 cm (4.9 GHz; one above bottom), and 20 cm (1.5 GHz; bottom). 



